Back

Chemical Science

Royal Society of Chemistry (RSC)

All preprints, ranked by how well they match Chemical Science's content profile, based on 71 papers previously published here. The average preprint has a 0.07% match score for this journal, so anything above that is already an above-average fit. Older preprints may already have been published elsewhere.

1
Fibril Paint: a class of amyloid-targeting peptides

Pedrola, J. A.; Dekker, F. A.; Guttmann, K.; van Leeuwen, L. M.; Singh, S.; Mayer, G.; Garfagnini, T.; Friedler, A.; Rudiger, S. G. D.

2024-08-25 biochemistry 10.1101/2024.08.25.609586 medRxiv
Top 0.1%
22.9%
Show abstract

Amyloid fibrils are a common pathological hallmark in multiple neurodegenerative diseases, yet molecular tools to selectively recognise and manipulate them remain scarce. We report FibrilPaints, a family of modular peptides designed for selective amyloid binding and adaptable chemical functionality. The degenerative amyloid-targeting unit of FibrilPaints, W5P4H3R2, has a high content of {pi}-stacking and aromatic side chains. Systematic sequence variation, altering charge, termini, and residue order, revealed the importance of the composition of the amyloid-targeting unit for high-affinity binding across Tau and Huntingtin fibrils. Importantly, sequence changes outside this unit do not preclude fibril binding, which permits attachment of fluorophores or E3-recruiting motifs for targeted protein degradation. This work establishes FibrilPaint as a modular peptide system for the detection and modulation of amyloids. O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=84 SRC="FIGDIR/small/609586v4_ufig1.gif" ALT="Figure 1"> View larger version (21K): org.highwire.dtl.DTLVardef@1e9cadborg.highwire.dtl.DTLVardef@50f9eorg.highwire.dtl.DTLVardef@446ee8org.highwire.dtl.DTLVardef@cd47a2_HPS_FORMAT_FIGEXP M_FIG Figure of content The modular design of FibrilPaints enables systematic evaluation of their functionality by testing the binding capacity of each variant (FibrilPaintX) to distinct amyloid fibrils. Successful binding results in visible painting of the fibrils, facilitating their detection and downstream research. C_FIG

2
Structure and dynamics of Burkholderia pseudomallei OXA-57, a distinctive low-efficiency class D β-lactamase with carbapenem-hydrolyzing activity

Bragginton, E.; Colenso, C.; Calvopina, K.; Hinchliffe, P.; Shaw, J.; Tooke, C.; Seng, R.; Chantratita, N.; Mulholland, A. J.; Schofield, C.; Spencer, J.

2024-12-23 biochemistry 10.1101/2024.12.23.630153 medRxiv
Top 0.1%
22.6%
Show abstract

The Gram-negative bacterium Burkholderia pseudomallei causes the severe disease melioidosis. {beta}-Lactams, including carbapenems, are the primary treatment, but are susceptible to chromosomal {beta}-lactamases, including the class D enzyme OXA-57. Here we show recombinant OXA-57 is active towards penicillins and first-generation cephalosporins, slowly hydrolyzes carbapenems including imipenem and meropenem, but is inactive towards oxyimino-cephalosporins (e.g., ceftazidime). Unlike many OXA enzymes, OXA-57 is sensitive to the mechanism-based inhibitor clavulanic acid, but less so to the diazabicyclooctane avibactam and not to the cyclic boronate vaborbactam. Crystal structures of covalent OXA-57:avibactam and OXA-57:meropenem complexes reveal limited hydrogen- bonding interactions with bound ligands. In molecular dynamics simulations, bound meropenem is mobile, while the water necessary for deacylation has only limited active site access. These observations are consistent with the low level of meropenem turnover, supporting proposals that OXA {beta}-lactamases generally possess limited carbapenemase activity, and highlight the potential importance of OXA-57 in B. pseudomallei {beta}-lactam resistance. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=109 SRC="FIGDIR/small/630153v1_ufig1.gif" ALT="Figure 1000"> View larger version (24K): org.highwire.dtl.DTLVardef@d635d7org.highwire.dtl.DTLVardef@168d67borg.highwire.dtl.DTLVardef@1d95733org.highwire.dtl.DTLVardef@1ebeabb_HPS_FORMAT_FIGEXP M_FIG C_FIG

3
In Silico Design And Binding Mechanism Of E3 Ligase Ubr1 Recruiters

Maria Solano, M. A.; Lazim, R.; Choi, S.

2024-01-15 biophysics 10.1101/2024.01.14.575594 medRxiv
Top 0.1%
19.6%
Show abstract

Proteolysis Targeting Chimeric Molecules (PROTACs) represent a promising avenue in drug discovery, as they can induce the targeted degradation of disease-relevant proteins within the cellular machinery. These compounds comprise a ligand tailored to bind the specific targeted protein connected to a recruiter molecule that engages with the E3 ligase. Despite their promise as therapeutic agents, their clinical advancement has encountered substantial challenges, primarily due to the limited availability of suitable E3 ligases. Additionally, cell permeability and proteolytic stability due to their peptide nature often hamper their application. In this study, we focus on the development of recruiters for the E3 ligase UBR1. This widely expressed protein has recently been demonstrated to be efficient in driving the degradation of oncogenic proteins. Our computational approach leverages a fragment-based peptidomimetics strategy, integrating pharmacophore filtering, docking, and fragment-linking optimization. Finally, we subject the wild-type peptide and the most promising combined fragments to advanced binding free energy calculations, unveiling insights into their dynamic water-mediated binding mechanisms and their potential as robust E3 ligase UBR1 recruiters, ultimately leading to the identification of promising compounds. This computational workflow is readily applicable to the development of related PROTACs and also to model protein-protein interactions with similar characteristics.

4
Molecular insights into inhibitor action on a key bacterial metabolic enzyme Cystathionine β-Synthase

Polepalli, S.; Roy, A.; Mondal, B.; Singh, A.; Dutta, S.

2025-07-05 biophysics 10.1101/2025.07.03.662948 medRxiv
Top 0.1%
19.5%
Show abstract

Tuberculosis (TB) remains a major global health threat, with Mycobacterium tuberculosis (Mtb) infecting nearly a quarter of the global population. Drug-resistant TB and HIV-TB co- infections emphasize the need for novel therapeutic approaches targeting essential metabolic pathways. Here, we investigated Mtb cystathionine {beta}-synthase (MtbCBS), a PLP-dependent enzyme critical for sulfur metabolism and redox regulation, owing its potential as a therapeutic target. We present the first high-resolution cryo-EM structure of MtbCBS bound to aminooxyacetic acid (AOAA), and employed molecular dynamics (MD) simulations, quantum mechanics/molecular mechanics (QM/MM) calculations, and comparative inhibition studies to reveal the molecular basis and determinants governing irreversible PLP-enzyme inhibition. Our Cryo-EM structural analysis revealed two highly conserved active-site residues, T75 and Q147, critically stabilizing the inhibitor complex. Through molecular mimic studies, we demonstrated the precise structural factors and electronic features critical for inhibition efficiency. These findings provide the first mechanistic rationale for PLP enzyme inhibition and offer a generalizable framework for designing covalent inhibitors targeting PLP-dependent enzymes implicated in infectious diseases, cancer, and neurological disorders. O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=115 SRC="FIGDIR/small/662948v1_ufig1.gif" ALT="Figure 1"> View larger version (35K): org.highwire.dtl.DTLVardef@195c7c7org.highwire.dtl.DTLVardef@6ee0c8org.highwire.dtl.DTLVardef@1679dfborg.highwire.dtl.DTLVardef@1b4e000_HPS_FORMAT_FIGEXP M_FIG O_FLOATNOGraphical Abstract.C_FLOATNO Proposed mechanism of MtbCBS inhibition by AOAA. Proposed mechanism of action of AOAA illustrated through a graphical representation. The cartoon shows the domains of MtbCBS. Tetrameric structure of MtbCBS with active site highlighted (in red dashes). Enlarged view indicating the PLP-AOAA blocked intermediate interacting with residue Q147. Reaction scheme depicting the previous missing understanding of active site residue role in -elimination. Schematic representation of AOAA inhibition mechanism in the catalytic cycle of MtbCBS. Reactions were drawn using ChemDraw. The proposed scheme highlights that AOAA (highlighted Red) binding is stabilized by nearby residues Q147 and T75, leading to the formation of PLP-AOAA covalent adduct, and substrates (highlighted green) were incapable of reverting the suicidal inhibition. The directionality and reversibility were mentioned for each step. C_FIG

5
Sequence-dependent transferability of the LRLLR membrane translocation motif: A computational study of smacN and NR2B9c peptides.

Munoz-Gacitua, D.; Blamey, J.

2026-02-19 bioinformatics 10.64898/2026.02.19.706823 medRxiv
Top 0.1%
19.0%
Show abstract

The LRLLR cell-penetrating motif can be transferred to confer membrane translocation activity, but only to compatible recipient peptides. Using umbrella sampling molecular dynamics simulations, we demonstrate that C-terminal LRLLR addition to the pro-apoptotic smacN peptide eliminates its translocation barrier entirely, transforming a +65 kJ/mol barrier into a -50 kJ/mol energy well. In contrast, N-terminal LRLLR addition to the neuroprotective NR2B9c peptide increases the translocation barrier from +85 to +100 kJ/mol, demonstrating that motif transfer can prove counterproductive for incompatible sequences. Cell-penetrating peptides offer promising strategies for intracellular delivery of therapeutic cargo, yet the sequence determinants governing their activity remain incompletely understood. The LRLLR motif, identified through systematic screening as essential for spontaneous membrane translocation, represents a minimal penetrating element whose transferability has not been previously evaluated. We appended this motif to two clinically relevant peptides: smacN, a tetrapeptide targeting inhibitor of apoptosis proteins in chemotherapy-resistant cancers, and NR2B9c, a nonapeptide that disrupts excitotoxic signaling in ischemic stroke. Potential of mean force profiles calculated across a POPC/POPG bilayer, combined with analysis of hydrogen bonding patterns, secondary structure propensity, and conformational dynamics, reveal the structural basis for these divergent outcomes. Successful transfer to smacN results from favorable complementarity: the hydrophobic, neutral smacN provides an ideal platform for the charged, amphipathic LRLLR motif, yielding a chimera capable of simultaneous interaction with both membrane leaflets. Transfer failure with NR2B9c stems from conformational rigidity induced by intramolecular hydrogen bonding, which prevents optimal membrane insertion, combined with unfavorable positioning of internal polar residues at the bilayer center. These findings establish that cell-penetrating motif transfer requires compatibility in charge distribution, hydrophobicity, and conformational flexibility between the motif and recipient sequence. The smacN-LRLLR chimera emerges as a promising candidate for experimental validation as a membrane-permeable therapeutic for survivin-positive tumors. More broadly, this work demonstrates the value of computational screening to identify compatible motif-cargo pairings prior to experimental investment.

6
High Throughput Virtual Screening and Validation of a SARS-CoV-2 Main Protease Non-Covalent Inhibitor

Clyde, A.; Galanie, S.; Kneller, D. W.; Ma, H.; Babuji, Y.; Blaiszik, B.; Brace, A.; Brettin, T.; Chard, K.; Chard, R.; Coates, L.; Foster, I.; Hauner, D.; Kertesz, V.; Kumar, N.; Lee, H.; Li, Z.; Merzky, A.; Schmidt, J. G.; Tan, L.; Titov, M.; Trifan, A.; Turilli, M.; Van Dam, H.; Chennubhotla, S. C.; Jha, S.; Kovalevsky, A.; Ramanathan, A.; Head, M.; Stevens, R.

2021-03-27 biophysics 10.1101/2021.03.27.437323 medRxiv
Top 0.1%
18.5%
Show abstract

Despite the recent availability of vaccines against the acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the search for inhibitory therapeutic agents has assumed importance especially in the context of emerging new viral variants. In this paper, we describe the discovery of a novel non-covalent small-molecule inhibitor, MCULE-5948770040, that binds to and inhibits the SARS-Cov-2 main protease (Mpro) by employing a scalable high throughput virtual screening (HTVS) framework and a targeted compound library of over 6.5 million molecules that could be readily ordered and purchased. Our HTVS framework leverages the U.S. supercomputing infrastructure achieving nearly 91% resource utilization and nearly 126 million docking calculations per hour. Downstream biochemical assays validate this Mpro inhibitor with an inhibition constant (Ki) of 2.9 {micro}M [95% CI 2.2, 4.0]. Further, using room-temperature X-ray crystallography, we show that MCULE-5948770040 binds to a cleft in the primary binding site of Mpro forming stable hydrogen bond and hydrophobic interactions. We then used multiple {micro}s-timescale molecular dynamics (MD) simulations, and machine learning (ML) techniques to elucidate how the bound ligand alters the conformational states accessed by Mpro, involving motions both proximal and distal to the binding site. Together, our results demonstrate how MCULE-5948770040 inhibits Mpro and offers a springboard for further therapeutic design. O_TEXTBOXSignificance StatementThe ongoing novel coronavirus pandemic (COVID-19) has prompted a global race towards finding effective therapeutics that can target the various viral proteins. Despite many virtual screening campaigns in development, the discovery of validated inhibitors for SARS-CoV-2 protein targets has been limited. We discover a novel inhibitor against the SARS-CoV-2 main protease. Our integrated platform applies downstream biochemical assays, X-ray crystallography, and atomistic simulations to obtain a comprehensive characterization of its inhibitory mechanism. Inhibiting Mpro can lead to significant biomedical advances in targeting SARS-CoV-2 treatment, as it plays a crucial role in viral replication. C_TEXTBOX

7
Time-resolved ligand dynamics revealed in a β-lactamase using room-temperature serial crystallography

Hinchliffe, P.; Tooke, C. L.; Beer, M.; Aller, P.; Kamps, J. J.; Parkinson, L.; Zhou, T.; Devenish, N.; Gu, D.-H.; Shilova, A.; Rodrigues, M. J.; Telfer, A.; Butryn, A.; Hanzevacki, M.; Freeman, E. I.; Park, S.; Schofield, C. J.; Park, J.; Owen, R. L.; Mulholland, A. J.; Orville, A. M.; Spencer, J.

2025-12-19 biochemistry 10.64898/2025.12.17.692331 medRxiv
Top 0.1%
18.5%
Show abstract

{beta}-Lactamases catalyze {beta}-lactam antibiotic hydrolysis and are important contributors to bacterial antimicrobial resistance; {beta}-lactamase inhibitors are widely used to overcome {beta}-lactamase-mediated antibiotic resistance. Nucleophilic serine {beta}-lactamases (SBLs) react with their substrates and clinically available inhibitors via a covalent reaction to give complexes which can undergo further reaction. Using room temperature drop on fixed target serial crystallography, where ligands are rapidly mixed with microcrystals, and classical single-crystal crystallography at cryogenic temperatures, we investigate the reversible covalent reaction of the SBL CTX-M-15 with the diazobicyclooctane inhibitor avibactam. We observe avibactam covalently reacted (ring-opened) with the nucleophilic Ser70, at timepoints from 80 ms to minutes (room temperature) and hours (100 K). These crystallographic data reveal time-dependent movement of the avibactam carbamoyl complex, from 1.3 s onwards, that has implications for the 5-exo-trig recyclization mechanism that determines inhibitor reformation. Combined with molecular dynamics simulations and quantum mechanics calculations at the density functional theory level, the results show that in the first seconds of the reaction the avibactam N-sulfate nitrogen is poorly positioned for recyclization. This subsequently equilibrates after 10 s to a stable endpoint that is in a conformation potentially primed to initiate recyclization through attack of the N-sulfate nitrogen on the carbamoyl carbon. These results further demonstrate the capacity of room-temperature serial crystallography to capture time-resolved changes in ligand conformation at an enzyme active site, complementing discrete classical cryo-crystallography. These data inform on ligand dynamics and the stereoelectronics of diazobicyclooctane inhibition, aiding drug discovery efforts to develop inhibitors of nucleophilic enzymes.

8
New from old: discovery of the novel antibiotic actinomycin L in Streptomyces sp. MBT27

Machushynets, N.; Elsayed, S. S.; Du, C.; Siegler, M. A.; de la Cruz, M.; Genilloud, O.; Hankemeier, T.; van Wezel, G. P.

2021-10-12 microbiology 10.1101/2021.10.12.464064 medRxiv
Top 0.1%
18.4%
Show abstract

Streptomycetes are major producers of bioactive natural products, including the majority of the antibiotics. While much if the low-hanging fruit has been discovered, it is predited that less than 5% of the chemical space has been mined. Here, we describe the novel actinomycins L1 and L2, which are produced by Streptomyces sp. MBT27. The molecules were discovered via metabolic analysis combined with molecular networking of cultures grown with different combinations of carbon sources. Actinomycins L1 and L2 are diastereoisomers, and the structure of actinomycin L2 was resolved using NMR and single crystal X-ray crystallography. Actinomycin L is formed via a unique spirolinkage of anthranilamide to the 4-oxoproline moiety of actinomycin X2, prior to the condensation of the actinomycin halves. Feeding anthranilamide to cultures of Streptomyces antibioticus, which has the same biosynthetic gene cluster as Streptomyces sp. MBT27 but only produces actinomycin X2, resulted in the production of actinomycin L. This shows that actinomycin L results from joining two distinct metabolic pathways, namely those for actinomycin X2 and for anthranilamide. Actinomycins L1 and L2 showed significant antimicrobial activity against Gram- positive bacteria. Our work shows how new molecules can still be identified even in the oldest of natural product families. IMPORTANCEActinomycin was the first antibiotic discovered in an actinobacterium by Selman Waksman and colleagues, as early as 1940. This period essentially marks the start of the golden era of antibiotic discovery. Over time, emerging antimicrobial resistance (AMR) and the declining success rate of antibiotic discovery resulted in the current antibiotic crisis. We surprisingly discovered that under some growth conditions, Streptomyces sp. MBT27 can produce actinomycins that are significantly different from those that have been published so far. The impact of this work is not only that we have discovered a novel molecule with very interesting chemical modifications in one of the oldest antibiotics ever described, but also that this requires the combined action of primary and secondary metabolic pathways, namely the biosynthesis of anthranilamide and of actinomycin X2, respectively. The implication of the discovery is that even the most well-studied families of natural products may still have surprises in store for us.

9
Cyclic peptide structure prediction and design using AlphaFold

Rettie, S. A.; Campbell, K. V.; Bera, A. K.; Kang, A.; Koslov, S.; De La Cruz, J.; Adebomi, V.; Zhou, G.; DiMaio, F.; Ovchinnikov, S.; Bhardwaj, G.

2023-02-26 biochemistry 10.1101/2023.02.25.529956 medRxiv
Top 0.1%
18.4%
Show abstract

Deep learning networks offer considerable opportunities for accurate structure prediction and design of biomolecules. While cyclic peptides have gained significant traction as a therapeutic modality, developing deep learning methods for designing such peptides has been slow, mostly due to the small number of available structures for molecules in this size range. Here, we report approaches to modify the AlphaFold network for accurate structure prediction and design of cyclic peptides. Our results show this approach can accurately predict the structures of native cyclic peptides from a single sequence, with 36 out of 49 cases predicted with high confidence (pLDDT > 0.85) matching the native structure with root mean squared deviation (RMSD) less than 1.5 [A]. Further extending our approach, we describe computational methods for designing sequences of peptide backbones generated by other backbone sampling methods and for de novo design of new macrocyclic peptides. We extensively sampled the structural diversity of cyclic peptides between 7-13 amino acids, and identified around 10,000 unique design candidates predicted to fold into the designed structures with high confidence. X-ray crystal structures for seven sequences with diverse sizes and structures designed by our approach match very closely with the design models (root mean squared deviation < 1.0 [A]), highlighting the atomic level accuracy in our approach. The computational methods and scaffolds developed here provide the basis for custom-designing peptides for targeted therapeutic applications.

10
Probing the energy landscape of α-Synuclein amyloid fibril formation by systematic K-to-Q mutagenesis

Kunka, A.; Farzadfard, A.; Larsen, J. A.; Saraceno, F.; Norrild, R. K.; Fricke, C.; Mohammad-Beigi, H.; Sadek, A.; Folke, J.; Aznar, S.; Buell, A. K.

2025-11-03 biophysics 10.1101/2025.11.01.685997 medRxiv
Top 0.1%
18.4%
Show abstract

The aggregation of natively disordered -Synuclein (Syn) into amyloid fibrils is a hallmark of Parkinsons and other neurodegenerative diseases. Understanding Syns pathological role remains a major challenge due to its complex, context-dependent energy landscape characterized by conformational plasticity and fibril polymorphism. Here, we present a systematic mutational analysis as a quantitative probe of the Syn energy landscape, focusing on electrostatic contributions to key aggregation pathways. We engineered Syn variants with one to eight lysine-to-glutamine substitutions and analyzed their aggregation under controlled conditions to delineate their effects on nucleation, elongation, seed amplification, fibril stability, and fibril polymorphism. We find that Syn aggregation from a homogenous solution can be modelled well using global properties, including protein concentration, charge, and ionic strength. Microscopic pathways and the resulting fibril polymorphs are instead modulated by sequence-specific effects. We identify mutations of residues found in fibril cores as perturbations that significantly modify the Syn free energy landscape, creating pathways and energy minima not accessible to the WT under the same experimental conditions. In contrast, mutations outside of the fibril core affect the magnitude of the relevant energy barriers whilst overall maintaining a WT-like free energy landscape. Our work outlines a scalable, quantitative framework that increases the informational output of the mutational studies of Syn using conventional assays. The approach can be extended by incorporating additional mutational and functional data to deepen our understanding of Syns energy landscape and its role in health and disease.

11
Nanoliter-scale selection of optimized bioengineered peptide antibiotics that rescue mice with bacterial lung infection

Boehringer, N.; Wuisan, Z. G.; Marner, M.; Kresna, D.; Mettal, U.; Schmitt, S.; Reiter, S.; Liu, Y.; Brinkrolf, K.; Rupp, O.; Schwengers, O.; Findeisen, J. I.; Herold, S.; Matt, U.; Schaeberle, T. F.

2024-06-01 microbiology 10.1101/2024.05.30.596569 medRxiv
Top 0.1%
18.3%
Show abstract

Increasing numbers of multi-drug resistant pathogens call for new chemical scaffolds, addressing novel targets, that can serve as lead structures for the development of life-saving drugs. For antibiotics, natural product-inspired molecules represent a most promising resource. Natural products evolved to high chemical complexity and occupy a chemical space different than synthetic libraries. However, clinical translation of promising natural products is often impeded by their relative inaccessibility to medicinal chemistry optimization, e.g. iterative synthesis of large series of derivatives. Here, this limitation is addressed with a randomized library of bicyclic heptapeptides based on the natural product darobactin that hits the clinically not addressed target BamA. Variants of the ribosomally synthesized and post-translationally modified peptides were generated using heterologous mutasynthesis. A parallelized screening assay is adapted in nanoliter-scale beads to test the darobactin derivatives against our sensor strain. Loss of fluorescence sorting prioritized 563 events out of the analyzed [~]500k beads. Re-testing confirmed 48 hit events, of which 40 proved to produce distinct darobactin-type molecules. Most promising structures were isolated and the growth inhibitory effects against Gram-negative pathogens validated. One of our current frontrunner compounds (i.e., darobactin B) was reinforced by the randomized screen. While microbiological investigations of the new derivatives is ongoing, darobactin B was profiled in later tier assays and compared to another promising, rationally-designed analog (i.e., darobactin B9, "D22"). Early ADMET profiling and efficacy tests in a mouse pneumonia model were performed. Darobactin B reduced bacterial load of Pseudomonas aeruginosa and Klebsiella pneumoniae by intraperitoneal, as well as intratracheal administration. Our study showcases the potential of mutasynthetic libraries for high-throughput screening and identification of functional peptides for drug lead discovery.

12
Molecular dynamics simulations demonstrate reduced antibiotic affinity to mirror bacterial targets

Fady, P.-E.; Ciccone, J.

2025-12-17 biochemistry 10.64898/2025.12.17.694932 medRxiv
Top 0.1%
18.3%
Show abstract

"Mirror life", self-replicating organisms composed of non-natural-chirality biomacromolecules, presents a severe future threat with potentially global consequences. Consequently, there is strong agreement among experts that it should not be created. However, there is some disagreement over how effective existing medical countermeasures might prove against mirror bacteria in the event that they were created. Here, we leverage computational chemistry methods including docking and molecular dynamics to determine the likely binding efficacy of existing antibiotics against natural and mirror bacterial protein targets. We find that existing most antibiotics fail to bind to mirror bacterial protein targets, unlike their natural chirality targets. This suggests that current medical countermeasures would not successfully exert an antimicrobial activity against mirror bacteria if the latter were created. Our results motivate further policy advocacy to curtail research that directly leads to the creation of mirror life.

13
Mechanistic Insight from Full Quantum Mechanical Modeling: Laccase as a Detoxifier of Aflatoxins

Zaccaria, M.; Dawson, W.; Kish, D. R.; Reverberi, M.; Bonaccorsi di Patti, M. C.; Domin, M.; Cristiglio, V.; Dellafiora, L.; Gabel, F.; Nakajima, T.; Genovese, L.; Momeni, B.

2021-01-03 biophysics 10.1101/2020.12.31.424992 medRxiv
Top 0.1%
18.3%
Show abstract

We demonstrate a path towards full Quantum Mechanics (QM) characterization of enzymatic activity. As a case-study, we investigate the detoxification of aflatoxin, a carcinogenic food contaminant, by laccase, a versatile oxidase capable of--but not efficient for--degrading aflatoxin. We use a combination of quantitative experimentation and QM modeling to show that low enzymatic steric affinity for aflatoxin is the main bottleneck, rather that the oxidative activity of laccase. To identify the structural elements responsible for low reaction rates, we perform a density functional theory (DFT) based modeling of both the substrate and the enzyme in a full QM simulation of more than 7,000 atoms. Thanks to our approach we point to amino acid residues that determine the affinity of laccase for aflatoxin. We show that these residues are substrate-dependent, making a full QM approach necessary for enzyme optimization. Altogether, we establish a roadmap for rational enzyme engineering applicable beyond our case study.

14
AMPGAN v2: Machine Learning Guided Discovery of Anti-Microbial Peptides

Van Oort, C. M.; Ferrell, J. B.; Remington, J. M.; Wshah, S.; Li, J.

2020-11-20 biophysics 10.1101/2020.11.18.388843 medRxiv
Top 0.1%
18.1%
Show abstract

Antibiotic resistance is a critical public health problem. Each year ~2.8 million resistant infections lead to more than 35,000 deaths in the U.S. alone. Antimicrobial peptides (AMPs) show promise in treating resistant infections. However, applications of known AMPs have encountered issues in development, production, and shelf-life. To drive the development of AMP-based treatments it is necessary to create design approaches with higher precision and selectivity towards resistant targets. Previously we developed AMPGAN and obtained proof-of-concept evidence for the generative approach to design AMPs with experimental validation. Building on the success of AMPGAN, we present AMPGAN v2 a bidirectional conditional generative adversarial network (BiCGAN) based approach for rational AMP design. AMPGAN v2 uses generator-discriminator dynamics to learn data driven priors and controls generation using conditioning variables. The bidirectional component, implemented using a learned encoder to map data samples into the latent space of the generator, aids iterative manipulation of candidate peptides. These elements allow AMPGAN v2 to generate of candidates that are novel, diverse, and tailored for specific applications--making it an efficient AMP design tool.

15
In Silico Study of the Early Stages of Aggregation of β-Sheet Forming Antimicrobial Peptide GL13K

Niknam Hamidabad, M.; Watson, N. A.; Wright, L. N.; Mansbach, R. A.

2024-01-28 biophysics 10.1101/2024.01.25.577308 medRxiv
Top 0.1%
18.0%
Show abstract

Antimicrobial peptides (AMPs) are of growing interest as potential candidates for antibiotics to which antimicrobial resistance increases slowly. In this article, we perform the first in silico study of the synthetic {beta} sheet-forming AMP GL13K. Through atomistic simulations of single and multipeptide systems under different charge conditions, we are able to shine a light on the short timescales of early aggregation. We find that isolated peptide conformations are primarily dictated by sequence rather than charge, whereas changing charge has a significant impact on the conformational free energy landscape of multi-peptide systems. We demonstrate that the lack of charge-charge repulsion is a sufficient minimal model for experimentally observed aggregation. Overall, our work explores the molecular biophysical underpinnings of the first stages of aggregation of a unique AMP, laying necessary groundwork for its further development as an antibiotic candidate.

16
De novo design of parallel and antiparallel A3B3 heterohexameric alpha-helical barrels

Chubb, J. J.; Albanese, K. I.; Rodger, A.; Woolfson, D. N.

2024-09-25 synthetic biology 10.1101/2024.09.25.614870 medRxiv
Top 0.1%
17.8%
Show abstract

The de novo design of -helical coiled-coil peptides is advanced. Using established sequence-to-structure relationships, it is possible to generate various coiled-coil assemblies with predictable numbers and orientations of helices. Here we target new assemblies, namely A3B3 heterohexamer -helical barrels. These designs are based on pairs of sequences with 3-heptad repeats (abcdefg) programmed with a = Leu, d = Ile, e = Ala, and g = Ser, and b = c = Glu to make the acidic (A) chains and b = c = Lys in the basic (B) chains. These design rules ensure that the desired oligomeric state and stoichiometry are readily achieved. However, controlling the orientation of neighboring helices (parallel or anti-parallel) is less straightforward. Surprisingly, we find that assembly and helix orientation are sensitive to the starting position of the heptad repeats (the register) in the peptide sequences. Peptides starting at g (g-register) form a parallel 6-helix barrel in solution and in an X-ray crystal structure, whereas the b- and c-register peptides form an antiparallel complex. In lieu of experimental X-ray structures for b- and c-register peptides, AlphaFold-Multimer is used to predict atomistic models. However, considerably more sampling than the default value is required to match the models and the experimental data, as many confidently predicted and plausible models are generated with incorrect helix orientations. This work reveals the previously unknown influence of heptad register on the orientation of -helical coiled-coil peptides and provides insights for the modeling of oligopeptide coiled-coil complexes with AlphaFold.

17
Tying the Knot: In Silico Design of Foldable Lasso Peptides

Nguyen, J. D. M.; da Hora, G. C. A.; Mifflin, M. C.; Roberts, A. G.; Swanson, J. M. J.

2025-01-22 biophysics 10.1101/2025.01.17.633674 medRxiv
Top 0.1%
17.6%
Show abstract

Lasso peptides are a unique class of natural products with distinctively threaded structures, conferring exceptional stability against thermal and proteolytic degradation. Despite their promising biotechnological and pharmaceutical applications, reported attempts to prepare them by chemical synthesis result in forming the nonthreaded branched-cyclic isomer, rather than the desired lassoed structure. This is likely due to the entropic challenge of folding a short, threaded motif prior to chemically mediated cyclization. Accordingly, this study aims to better understand and enhance the relative stability of pre-lasso conformations--the essential precursor to lasso peptide formation--through sequence optimization, chemical modification, and disulfide incorporation. Using Rosetta fixed backbone design, optimal sequences for several class II lasso peptides are identified. Enhanced sampling with well-tempered metadynamics confirmed that designed sequences derived from the lasso structures of rubrivinodin and microcin J25 exhibit a notable improvement in pre-lasso stability relative to the competing nonthreaded conformations. Chemical modifications to the isopeptide bond-forming residues of microcin J25 further increase the probability of pre-lasso formation, highlighting the beneficial role of non-canonical amino acid residues. Counterintuitively, the introduction of a disulfide cross-link decreased pre-lasso stability. Although cross-linking inherently constrains the peptide structure, decreasing the entropic dominance of unfolded phase space, it hinders the requisite wrapping of the N-terminal end around the tail to adopt the pre-lasso conformation. However, combining chemical modifications with the disulfide cross-link results in further pre-lasso stabilization, indicating that the ring modifications counteract the constraints and provide a cooperative benefit with cross-linking. These findings lay the groundwork for further design efforts to enable synthetic access to the lasso peptide scaffold. SIGNIFICANCELasso peptides are a unique class of ribosomally synthesized and post-translationally modified natural products with diverse biological activities and potential for therapeutic applications. Although direct synthesis would facilitate therapeutic design, it has not yet been possible to fold these short sequences to their threaded architecture without the help of biosynthetic enzyme stabilization. Our work explores strategies to enhance the stability of the pre-lasso structure, the essential precursor to de novo lasso peptide formation. We find that sequence design, incorporating non-canonical amino acid residues, and design-guided cross-linking can augment stability to increase the likelihood of lasso motif accessibility. This work presents several strategies for the continued design of foldable lasso peptides.

18
Acetylation of L-leucine switches its carrier from the L-amino acid transporter (LAT) to organic anion transporters (OAT)

Churchill, G. C.; Strupp, M.; Bremova-Ertl, T.; Factor, M.; Patterson, M. C.; Platt, F. M.; Galione, A.

2020-11-13 pharmacology and toxicology 10.1101/2020.11.12.379701 medRxiv
Top 0.1%
17.2%
Show abstract

N-acetyl-DL-leucine is an analogue of the alpha amino acid leucine with a chiral stereocenter. The active L-enantiomer of the racemate is currently under development for rare neurological disorders. Here we present evidence that a selective recognition of N-acetyl-L-leucine versus L-leucine by different uptake transporters significantly contributes to the therapeutic effects of N-acetyl-L-leucine. A previous study of the pharmacokinetics of racemic N-acetyl-DL-leucine and N-acetyl-L-leucine revealed D-L enantiomer competition and saturation kinetics, best explained by carrier-mediated uptake. The strategy we used was to first analyze the physicochemical properties associated with good oral bioavailable drugs and how these are alerted by N-acetylation by comparing N-acetyl-L-leucine with L-leucine. Using in silico computational chemistry we found that N-acetylation has a profound impact on certain physicochemical properties that can rationalize why N-acetyl-L-leucine is drug-like compared to L-leucine. Our calculations show that at physiological pH, L-leucine is a zwitterion, whereas N-acetyl-L-leucine is present as mainly an anion. Specifically, N-acetylation removes a charge from the nitrogen at physiological pH and N-acetyl-L-leucine is an anion that is then a substrate for the organic anion transporters. We examined N-acetyl-L-leucine uptake in human embryonic kidney cells overexpression candidate organic anion transporters (OAT) and pharmacological inhibitors. We found that N-acetyl-L-leucine is a translocated substrate for OAT1 and OAT3 with low affinity (Km ~10 mM). In contrast, L-leucine is known to be transported by the L-type Amino Acid Transporter (LAT) with high affinity (Km ~0.2 mM) and low capacity. The clinical consequence is that L-leucine uptake becomes saturated at 50-fold lower concentration than N-acetyl-L-leucine. These results demonstrate a mechanism of action that explains why N-acetyl-L-leucine is effective as a drug and L-leucine itself is not.

19
DEER of Singly Labelled Proteins to Evaluate Supramolecular Packing of Amyloid Fibrils

Tsay, K.; Equbal, A.; Li, Y.; Tsui, T.; Han, S.; Fichou, Y.

2026-02-09 biophysics 10.1101/2025.11.10.687719 medRxiv
Top 0.1%
17.1%
Show abstract

1The formation of protein amyloid fibrils in the brain is a hallmark of various neurodegenerative diseases, including Alzheimer diseases and Parkinson diseases. Amyloid fibrils are highly ordered aggregates in which proteins folded in two-dimensional layers stack along one axis to form elongated linear assemblies. The specific conformations adopted by the proteins within each layer of the amyloid correlates with the pathology. Furthermore, their spreading and templating competency relies on strict in register packing of the folded proteins along the fibril growing axis. There is a need for tools to characterize not only the protein fold across the fibril cross section but also the spatial ordering of the proteins stacked along the amyloid fibril axis. We present an approach based on double electron electron resonance spectroscopy (DEER) using singly labelled tau protein assembled in amyloid fibrils that can deliver an apparent dimensionality of the supramolecular organization of tau fibrils. The parameters of the DEER background function can be used to assess the amyloid core location and packing order, and track time-resolved formation of aggregation intermediates. Showcasing the method on tau, we demonstrate that heparin-induced tau fibrils are mispacked while seeded aggregation can template amyloid fibrils with a higher packing order. This study benchmarks a new method that will provide critical structural insights into amyloid assemblies. O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=134 SRC="FIGDIR/small/687719v2_ufig1.gif" ALT="Figure 1"> View larger version (36K): org.highwire.dtl.DTLVardef@8c42c9org.highwire.dtl.DTLVardef@74b0f8org.highwire.dtl.DTLVardef@10e8095org.highwire.dtl.DTLVardef@11bc727_HPS_FORMAT_FIGEXP M_FIG DEER-derived dimensionality parameter differentiates between spin-labelled tau protein in solution, in perfectly aligned fibril or in imperfectly aligned fibril (fibril structure edited from PDB 6qjh). C_FIG

20
Allostery links hACE2 binding, pan-variant neutralisation and helical extension in the SARS-CoV-2 Spike protein

Colyer, A.; Wolf, E.; Lento, C.; Ustav, M.; Sljoka, A.; Wilson, D. J.

2025-01-29 biochemistry 10.1101/2025.01.28.635280 medRxiv
Top 0.1%
17.1%
Show abstract

The SARS-CoV-2 spike protein is highly antigenic, with epitopes in three distinct regions of the receptor binding domain (RBD) alone that have known mechanisms of neutralization. In previous work, we predicted a fourth RBD epitope based on allosteric conformational perturbations measured by hydrogen-deuterium exchange mass spectrometry (HDX-MS) upon complexation with the canonical spike protein target, human angiotensin-converting enzyme 2 (hACE2). We subsequently identified a pan-neutralizing antibody (ICO-hu104) with the predicted epitope, however, as the epitope was somewhat distant from the hACE2 binding interface, and our previous work limited to the spike RBD, the neutralization mechanism was unclear. Using HDX-MS, we investigated the binding of ICO-hu104 to the full-length SARS-CoV-2 spike protein from Wuhan, Delta and Omicron variants. We demonstrate that binding of ICO-hu104 at its epitope results in an increase in deuterium uptake in the distant HR1 domain for variants of concern, which in a biological context could be indicative of destabilisation of the helices within this region, promoting S1 shedding or failure of helical extension during S2-mediated fusion. This is supported by our computational modelling, highlighting propagation of allosteric effects to the S2 coiled-coil region. Collectively, this work demonstrates an alternative neutralization mechanism for ICO-hu104 which is distinct from its first-generation predecessors and thus opens alternative avenues targeting non-RBD epitopes through assessment of allosteric perturbations. HighlightsO_LIHDX-MS reveals decreased deuterium uptake within the HR1 region of S2 for SARS-CoV-2 spike protein for variants of concern when bound to ICO-hu104. C_LIO_LIComputational modelling validates high allosteric coupling between ICO-hu104 epitope and HR1 region. C_LIO_LISuggests an alternative neutralization mechanism from its predecessor ICO-hu23, whereby destabilisation of helices within the HR1 region promotes S1 shedding and/or failure of helical extension during S2-mediated fusion. C_LI Graphical abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=81 SRC="FIGDIR/small/635280v1_ufig1.gif" ALT="Figure 1"> View larger version (20K): org.highwire.dtl.DTLVardef@1175bdborg.highwire.dtl.DTLVardef@8fb17borg.highwire.dtl.DTLVardef@1cd13f7org.highwire.dtl.DTLVardef@d9c121_HPS_FORMAT_FIGEXP M_FIG C_FIG